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The earlier-described cluster methods are used to generalize the Ginzburg-Landau gradient expansion for 
the free energy of an inhomogeneous alloy to the case of not small values of order parameters and variations of 
composition. The results obtained reveale a number of important differences with the expressions used in the 
phenomenological phase-field approach. Differential equations relating the local values of concentration and 
order parameters within the antiphase or interphase boundary (APB or IPB) are derived. These equations 
are applied to study the segregation at APBs near the phase transition lines; the structure of APBs and IPBs 
near tricritical points; wetting APBs in phases with single and several order parameters; and also some effects 
of anisotropy of APBs under Llo and Ll2-type orderings. 
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Studies of inhomogeneous alloys attract interest 
.from both fundamental and applied points of view, in 
particular, in connection with microstructure evolution 
under phase transformations [1-12]. Typical inhomo- 
geneities in such problems are antiphase or interphase 
boundaries (APBs or IPBs) which separate differently 
ordered domains or different phases. Both experimental 
and theoretical studies show that in situations of practi- 
cal interest the APB or IPB width usually much exceeds 
the interatomic distance [5-12]. Therefore, Ginzburg- 
Landau (GL) type gradient expansions can be used to 
describe the free energy of such states even though the 
order parameters and concentration variations here are 
typically not small, contrary to assumptions of the stan- 
dard GL theory. Employing such generalized GL func- 
tionals (suggested first by Cahn and Hilliard [0) is now 
referred to as the phase-field method, and it is widely 
■used for most different systems, see e. g. [6-8]. How- 
ever, a number of simplifying asumptions are usually 
employed in this phenomenological approach, and so its 
relation to more consistent theoretical treatments re- 
mains unclear. Recently microscopical cluster methods 
have been developed for inhomogeneous alloys [9-13]. 
Below I use these methods to derive the GL functionals 
and then apply them for studies of APBs and IPBs. 

To be definite, I consider a binary alloy A C B!_ C at 
c < 1/2. Various distributions of atoms over lattice sites 
i are described by the mean occupations c, = (rij) where 
rii is unity when the site i is occupied by atom A and 
zero otherwise, while averaging is taken, generally, over 
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the space- and time-dependent distribution function |9J . 
The free energy F{ci} in the cluster description can be 
written as a series MM: 



F 



E/^El^+E^+'-E^-' 



(i) 



Here f l is the free energy per site i; F{ = T[cj In Cj + 
(1 — Cj) ln(l — Cj)] is the mixing entropy contribution; 
pi...k _ F l ( c . t , ..cj,) is the contribution of interactions 
within /-site cluster of sites i, . . . fc; and m is the max- 
imum cluster size considered. The simplest mean-field 
approximation (MFA) and the pair-cluster one (PCA) 
correspond to neglecting many-site contributions F m> i 
m 

(0), while in a more refined, tetrahedron cluster ap- 
proximation - TCA (that should be used, in particular, 
to adequately describe the LI2 and Llo-type orderings 
§) Eq. (0) includes also 4-site terms F% kl @. 

For the homogeneous ordered structure, the mean 
occupation Cj = c(rj) at site j with the lattice vector 
Tj can be written as a superposition of concentration 
waves with some superstructure vectors k s [8-12]: 



7y s exp(ik s rj) = ^77,, exp(ik p r :) ). (2) 



Here amplitudes ?7 S can be considered as order parame- 
ters; the last expression includes also term with r) p = c 
and k p = 0; and for simplicity both parameters rj s and 
factors exp(ik s ri) are supposed to be real which is the 
case, in particular, for the B2, Llo and Ll2-type order. 
For weakly inhomogeneous states, amplitudes rj p in (0) 
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are not constants but smooth functions of coordinates 
IV Thus functions f{cj} in (Q) can be expanded in 
powers of differences Scj = Y^p^Vp ex P('^ s p r j) where 



pVp- 



(3) 



Here Vri p = drj p /drf, V Q/3 77 p = d 2 r) p /drf dr\ , r jt = 
{vj — Vi), and the summation over repeated Cartesian 
indices a,/3 — 1,2,3 is implied. After substitution of 
these expressions into Eq. (|l|) one can proceed from 
summation over i to integration over continuous vari- 
able r = rj. Making also standard manipulations with 
part-by-part integration of terms with V Qj ary p Q one 
obtains for the GL functional: 



F = 



ilpq v c 



V P VpVq + f{Vp] ■ ( 4 ) 



Here v a is volume per atom; f{rj p } is function f l {cj} 
in (Q) averaged over all sublattices with Cj given by Eq. 
(|2|); and is given by the expression: 



y pq 
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Sy exp[i(k p rj - k 9 r,)] 



where 5^ = d 2 f % j 'dckdcj . Note that the last term of 
(|^) is nonzero only when all three sites i, j and k are dif- 
ferent, and so it is present only when the non-pairwisc 
contributions F^yJ? in ([j]) are taken into account, such 
as the TCA terms F^ kl |||. 

For the LI2 and Llo phases in FCC alloys, Eq. (||) 
includes waves with three vectors k s : ki = (100)27r/a, 
k2 = (010)27r/a, and k3 = (001)27r/a where a is the 
lattice constant [8-15]. The local order within APBs 
in these phases can be described by the distribution 
of amplitudes of these waves (?7i, r?2j %) of the type 
(£, 77, 77) and corresponds to the tetragonal symmetry 
§, |, 0. To illustrate the form of terms in (|), 
we present their MFA and PCA expressions for this 
type local order. Tensors g p P here can be described 
in terms of their "transverse" and "anisotropic" compo- 
nents, g^ q = gf q = gf q and g a pq = (g£ - gf q ). Using for 
simplicity the 2-neighbor-interaction model: v n> 2 = 0, 
one obtains in the MFA: 



9(C = -\ a2v 2] 9am = ± \ a2v ^ 9^1 = l a2 ( v i - 2 «2); 



9cc = -^ 2 {vi + v 2 ); g a cc = g^ cM = 0, (6) 



9a,cc = ^(-fl ± -xt~ <P2 d )\ 9(c = \a 2 Vi ; 
5cc.cc = \<?{vt ± 9cc = \a 2 (-Vl - xt + 

± 1 2/ ab r, +\ a 1 2 ab 

9 m = 2 a Wi _ 2 *2 ); 9 m = <Pi ; 
9rt,vc = J^i^i ~ 2< ^2 ); 9$c;,r,c = ±^a 2 ipi- (7) 

Here plus or minus in (±) corresponds to the first or the 
second pair of lower indices in the left-hand-side of Eqs. 

®)-(0); 

= -Tf n /R» ; ^ = \ (<p£ ± <^ d ) ; 

V>± = ~ ± ^) ; x± = 5 (vC ± vi? ) ; 

i# = [1 + 2/ n (ci + Cj - 2c lCj ) + / 2 (c, - c,) 2 ] 1 / 2 ; (8) 

/„ = exp(—v n /T) — 1 is the Mayer function; and in- 
dex i or j equal to a, 6 or d corresponds to the mean 
occupation c, or Cj of one of three different sublattices: 

c a = c + ( + 2 V ; c b = c + (-2ri; c d = c-(. (9) 

If the TCA is used, the nearest-neighbor contributions 
ipY in Eqs. (0) are replaced by the relevant tetrahedron 
contributions presented in Ref. For the B2 or- 
der, there is only one order parameter 7y s = 77 |9) ; terms 
g p ^ have a cubic symmetry: g p P = S a /3g pq ; and the 
MFA and PCA expressions for g pq are similar to those 
for g± in Eqs. (§ and (@). 

In the phase field method, terms g p $ are assumed 
to not depend on local parameters r) r , and so to be zero 
at p 7^ q "on considerations of symmetry" [6-8]. Eqs. 
(||)-(0) show that it may correspond only to the simplest 
MFA, while in more accurate approaches, such as PCA 
and TCA, the dependences g pq (i]r) can be significant. 

Let us now consider the case of a plane APB (or 
IPB) when parameters r] p in (g) depend only on the dis- 
tance £ = rn where n = (cos a, sin a cosy), sin a sin <p) 
is normal to the APB plane. To find the equilibrium 
structure, one should minimize the functional (Q) over 
functions r) p (£) at the fixed total number of atoms ||. 
Let us first consider the APB between two B2-ordered 
domains. Then variational equations for the order pa- 
rameter 7y(£) and the concentration c(£) have the form: 

1 



g v W+9r > cc" + -g^(v') 2 +9T , v'c + 



(92 



cc\/„/\2 



9 c v c )(c 



while the PCA expressions for g^ and g p are: 



g V cv" + 9ccc" + « - -gTM) z + 

+ k cc = ^(/c-M)- 



(10) 



Ginzburg-Landau functionals for alloys with application to antiphase boundaries 



3 



Here prime means taking derivative over £; the lower 
index 77 or c means taking derivative over 77 or c; and /-i 
is the chemical potential. At £ — > oo functions c and 77 
tend to their equilibrium values, Co and r?o( c o)- 

Multiplying first and second Eq. (|l^) by 77' and c', 
respectively, summing them, and integrating the result, 
one obtains the first integral of this system of equations: 

9 m {il') 2 + 2 5i?c t/c' + g cc (c') 2 = Q(r], c) (11) 

where is the non-gradien part of the local excess grand 
canonical potential per atom: 



(12) 



and index zero at the function means its value at rj = r/ 
and c = Co . 

For what follows it is convenient to consider the or- 
der parameter 77 as an independent variable, while c and 
fl, as its functions determined by Eqs. (^0|)- (|l2|) . Then 
the dependence 77' (rj) is determined by Eq. (fTO): 



7/ = dr)/d£ = (Vl/G) 1/2 



(13) 



where G is (g vr , + 2g nc c + g cc c 2 ), and c is dc/dn. Using 
Eq. (|13|) one can exclude 77' from the system of equa- 
tions (|10| ) and obtain the differential equation for c(rj) 
to be called the composition-order equation (COE): 

[c(g 2 c -g cc g vv ) + ^]2n/G = 

= - fc) (9w + 9r,cC) + f n (g V c + 9ccC) (14) 

where $ is a linear function of derivatives g^ q : 



2 % 



$ = (.9»?c + 9ccc) 

~(9nv+ 9vcc){g7, c - 5 9, 



g™6- 



'/'/ 



9^+\gfc? 



(15) 



Because of the equilibrium conditions: /? = 0, f® = fi, 
function £l(r]) (|l2| ) at 77 — > 770 is proportional to (77— 770) 2 ■ 
Therefore, the initial value 6(770) can be found by taking 
the ?/ -> 770 limit of Eq. (Q). 

For the given solution 0(77) of COE, the coordinate 
dependence r)(£) is determined by integrating Eq. (|i"3|): 



e=a+ / ^(c/n) 1 / 2 



(16) 



where the reference point £1 is determined by the choice 
of value 771 = 77(^1). For the symmetrical APB for which 
77 — > ±770 at £ — > ±00, functions c, O and G are even 
in 77, and it is natural to put £1 = at 771 = 0. But 
for an IPB separating the ordered and the disordered 
phase, values 771 — > correspond to £ — > (—00), and so 
£1 should be chosen at some intermediate value 77!. 



The surface energy a and the surface segregation T 
is the excess of the grand canonical potential and of B 
atoms, respectively, per unit area 0). Taking into ac- 
count Eqs. one obtains for the surface energy: 



dr/ (Gfl) 1/2 



(17) 



where 77 m i n is (—770) for an APB and zero for an IPB, 
while the surface segregation is given by the expression: 



1 



d v (c -c) {G/n) 1 ' 2 . 



(18) 



For a symmetrical APB, the integral in (17) or ( |18| ) can 
be written as twice of the integral over positive 77. 

can also be derived for 



Relations similar to (11)- 
phases with several order parameters, such as the LI2 
or Llo phase. In particular, for an APB separating two 
Ll2-ordered domains, the order parameters (771,772,773) 
have the form (£, 77, 77) mentioned above with the lim- 
iting values (770,770,770) and (770,-770,-770) at £ — s- ±00. 
The variational equations and their first integral have 
the form analogous to Eqs. ( |l0| ) and ( |Tlj ) but include 
three functions, c(£), £(£) and T](^). It is again conve- 
nient to consider c and £ as functions of 77 (while for 
APBs in the Llo phase, c and 77 as functions of £ ) and 
obtain a system of equations for c(rf) and £(77) analo- 
gous to CO E ([14 ). Equations for £(77), a and T preserve 
their form (|l6|)-(i"§l) but G(rj) now includes the deriva- 
tive £ = d^/dr) and six functions g pq which are related 
to g^ a in Eqs. (§), (0) as follows: 



g pq (a) 



9pg + 9p 



cos 2 a, 



(19) 



where a is the angle between the APB orientation and 
the local tetragonality axis. When the nearest-neighbor 
interaction v\ much exceeds the rest ones (as in CuAu- 
based alloys [9-11]), the functions g pq (a) are highly 
anisotropic which is illustrated by Eqs. (^): g m ~ 
sin 2 a; g^ ~ cos 2 a. It results in a notable anisotropy 
of distributions of APBs, including the presence of many 
low-energy "conservative" APBs with a ~ in the LI2 
phase and a ~ 7r/2 in the Llo phase, as well as a pe- 
culiar alignment of APBs in "twinned" Llo structures 
[9-11]. In more detail, COE and Eqs. for the 

LI2 and L1q phases will be discussed elsewhere. 



Let us now discuss some applications of Eqs. (14)- 
(J]~g|) . First I consider the case when the equilibrium 
order parameter 770 is small. Then the function / in (^) 
can be written as the Landau expansion: 



7(77, c) = ip + arj 2 + br] + drf 



(20) 
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where ip, a, b and d are some functions of concentra- 
tion c and temperature T. The equilibrium value t/q is 
determined by the equation f® — 0, while the order- 
ing spinodal T = T s (c) (the disordered phase stability 
limit) is determined by the equation: a(c, T) = 0. Small 
values r/o under consideration correspond to cq, T points 
near the ordering spinodal where a(co,T) is small. 

It is clear from both physical considerations and the 
results below that the difference (c — Co) at small 770 is 
small, too. Therefore, functions (f c — /1) and f v in Eq. 
([14] ) can be expanded in powers of (c — cq), rj and i]o'- 



f c -H= a° c (v 2 - Vo) + Vcc (c - c ) + . 
fr, = 4?7 b {v 2 - Vo) + 2 V a° c (c - cq) 



(21) 
(22) 



Here ip°. c is (d 2 ip/dc 2 )o, and dots mean terms of higher 
orders in rj 2 ,. The analogous expansion of fi(c, 77) starts 
with terms bilinear in (c — Co) and (?7 2 — ?7q), while for 
functions g nc and c the expansions start with terms lin- 
ear in 77. Thus terms with f2 and / r; in (^) are propor- 
tional to ?7q being small compared to (f c — /j) ~ ?7q, and 
COE is reduced to the equation / c = /1, which yields: 



C0-cW-(% 2 -^)(-«c)/^c 



(23) 



Taking derivative of equation a(c,T) = one obtains: 
(-0°) = T' s a where a is {da/dT) and Tj = dT s /dc. 
Thus the surface segregation at APB is proportional to 
the ordering spinodal slope T s '(c ), and so it decreases 
with approaching the critical point where T' s = 0. 

Let us now suppose the alloy state Cq, T in the c, T 
plane to be close to the second-order transition line T s (c) 
far from possible tricritical points. Then the higher- 
order terms in expansions (p(i| )-(p2]) can be neglected, 
and for the function £l(r]) (|12[) such expansion yields: 



2\2. 



io-(«W. (24) 



Using Eqs. (|i~6|)-(|l8|) and ( p4] ) one obtains in this case 
for ?7(f;), c(£), the APB energy a and the segregation T: 

r)(Z) = 7jotanh(e/<5); c - c(£) = 770 Acosh- 2 ^); 
ff= 3^ %3 (^) 1/2; r = f *AGf/5)V» (25) 



1/2 



where A is aT'Jip° cc ; g is 5°^; and 5 = {g/r] 2 b) 
is the APB width. These expressions generalize the 
earlier MFA results |L2| to the case of any GL func- 
tional. The dependences 7?(£) and c(£) in Eqs. ( p5| ) 
are similar to those observed in the Monte Carlo study 
of segregation at APBs ||. The temperature depen- 
dence of this segregation at small rjo is more sharp than 
that of the APB energy: r oc rj ~ (T s - T) 1 / 2 while 
a oc T)l ~ (T s - Tf' 2 . Eqs. (§4|)-(§|) also show that 



the presence of segregation results in a renormalization 
of the Landau parameter & entering characteristics of 
APB to the lesser value b given by Eq. (^4|). It results 
in a decrease of the APB energy a and an increase of 
its width 6 and segregation T under decreasing temper- 
ature T along the ordering spinodal T — T s (c). 

The point co,T at which both a(co, T) in (|2p| ) and b 
in ( p4f ) vanish corresponds to the tricritical point Ct,T t . 
At T < T t the second-order transition line T s (c) in the 
c,T plane splits into two binodals, Cbo(T) and Cbd(T), 
delimiting the single-phase ordered and disordered field, 
respectively. Such tricritical point is observed, for exam- 
ple, in alloys Fe-Al At this point the lowest order 
terms in Eqs. (E^) and (j|J) vanish, one should con- 
sider the next-order terms, and function Q ([l2]) at small 
x = (cq — Cj) and t = (T — T t ) takes the form: 



ft( V )=A( V 2 -r 1 2 ) 2 (r, 2 + h). 



(26) 



Here h = h{x, t) is a linear function of x and t which can 
be written in terms of the binodal temperature deriva- 
tive c' bo = dcbo/dT as: h = v (x — tc' bo ), while A and v 
are some positive constants. Using Eqs. (|l~6|)-(|l8|) and 
(|26| ) one obtains for the characteristics of APB near T*: 

770 sinhy Ar/ 2 (l + a) 

f(0 = 7 — n — ; — znzi c o - c (0 - 



(cosh y + a) 1 / 2 



(cosh y + a)' 



J (a) 



4U/2. 



2XL(a) 



1/2 



(27) 



Here y is £/<5; a is rfc/h; while 8, L(a) and J(a) are: 

I - [g/A V 2 (h + V 2 o)] 1/2 ;L(a) = In [(1 + a) 1 ' 2 + a 1 ' 2 ]; 
J{a) = J-[(l + 4a)i(a) + (2a - l)(a + a 2 ) 1 / 2 ]. (28) 

Function h{x,t) in (p6|)-(p8|) is proportional to the dis- 
tance in the c, T plane from point co,T to the binodal 
Cbo{T), while i] 2 is proportional to the distance to the or- 
dering spinodal T s (c). Thus at small a«l, Eqs. (26)- 
dH) turn into with b = Ah and describe the 

critical behaviour of APB near T s (c) discussed above. 
The opposite case, a> 1, corresponds to the region of 
"wetting" APBs which received recently much attention 
H ||, |1| . Value h = corresponds to the ordered state 
with c = c bo (T), and then COE @ describes an IPB 
betwen this state and the disordered state with r/d = 
and Cd = Cbd(T). Substituting Eq. ( p6| ) with h = into 
Eqs. and (|l7|) one obtains for this IPB: 



no 



A% 2 



0"APB = 2<7ipb. (30) 
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Here z is f/ft; ft = {g / A) 1 ' 2 /2ril is the IPB width; 
ctipb = (.9-4??o) 1/2 /2w a is the IPB energy; c d is c - A^g; 
and the coordinate £1 = in ( |l6| ) is chosen at r\\ = 
f)o/v2 where 0(771) is (cq + c<j)/2. Eqs. (29) show that 
the order parameter 77 in the disordered phase decreases 
with moving from IPB much more slowly than the con- 
centration deviation: 77 ~ (c — c^) 1 / 2 . Eqs. (p7|)-(^0|) 
also show that in the "wetting" regime of large a, the 
profiles 77(5) and c(£) in Eq. ( p7|) correspond to the pres- 
ence at £± = ±ft In a of two almost independent IPBs 
described by Eqs. (29). The total width I ~ (f+ - f_) 
and the segregation F for such APB are proportional to 
Ln(l//i) while the energy difference (ctapb ~ 2ctipb) is 
proportional to h\n{l/h), which arc usual dependences 
for the wetting regime |l6). Eqs. (|27|)-(|30|) specify these 
relations for the vicinity of tricritical points and enable 
one to follow the transition from wetting to the critical 
behaviour of APBs under variation of T or cq. 

Let us now apply Eqs. (p^[)-(p7[) to derive the wet- 
ting relation ([}0|) for the phases with several order pa- 
rameters, in particular, for the Ll 2 phase in equilibrium 
with the Llo or the disordered FCC (Al) phase. This 
problem was discussed by a number of authors [|[ ^, |l5| 
but the general proof seems to be absent yet. Let us first 
note that in consideration of IPB, the initial condition 
to COE (|j^ ) can be put in either the ordered or the dis- 
ordered phase, i.e. at (c, 77) values equal to either (co, 770) 
or (cd,0), while the solution cipb(t7) at either choice is 
the same and unique. Therefore, in consideration of 
APB with the same initial values 00,770, the solution 
capb(?7) coincides with ctpb(t7) at 77 > 0, it is capb( — V) 
at 77 < 0, and so Eq. @) follows from Eq. (0). For the 
APB or IPB in the LI2 phase, the local order can be 
described by the parameters (c, £, 77) mentioned above, 
and their initial values in COE are (00,770,770), while 
the final ones are (00,770,-770), (c d ,0, 0), and (ci,rji,0) 
for the case of an APB, IPB(L1 2 -Al), and IPB(L1 2 - 
Llo), respectively, where q and 77; correspond to the 
second binodal LI2-LI0. Therefore, the wetting rela- 
tion j3(i| ) for the L1 2 -A1 or Ll 2 -Llo phase equilibrium 
follows from COE and Eq. ([17|), just as for the single- 
order-parameter case. Note, however, that at the given 
orientation no there are three types of an APB in the 
Ll 2 phase with the local order (771 , 772 , 773 ) of the form 
((■,77,77), (77, (,77) or (77,77,0, and the structure and en- 
ergy for each type is, generally, different. Therefore, 
there are at least three types of an IPB(L1 2 -A1) corre- 
sponding to "a half of the relevant APB in the wetting 
limit. In the course of the kinetical wetting (for ex- 
ample, under Al— sAl+Ll 2 transformations studied in 
J?], U) each APB first transforms into two "its own" 
IPBs, but later on these IPBs can evolve to other types. 



The effects of anisotropy under wetting APBs in Eqs. 
(|T(f)-(|Ts|) are described by a factor G 1 / 2 , while the main 
contribution to here is determined by the thermo- 
dynamic relations. In particular, singular contributions 
to a and T under wetting L1 2 -APB by the Al or Llo 
phase correspond to the region of small 77 where f2 has 
the same form as in Eq. (|26|), while G(rj) is reduced to 
its first term g vv as functions g^, g^c, ( and c, being 
odd in 77, vanish at small 77. Thus the main contributions 
to the APB width and energy take the form: 

I ~ g\j 2 In (l//i); ( ctapb - 2<np B ) ~ 9l'n Mn (V^O. 

where the angular dependence g m is given by (jl9|). 
Therefore, the wetting effects can reveale a significant 
anisotropy, particularly in the short-range interaction 
systems. It agrees with some previous results 1, |, ||. 

The author is much indebted to I.R. Pankratov for 
the help in this work. The work was supported by the 
Russian Fund of Basic Research under Grants No. 00- 
02-17692 and 00-15-96709. 
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